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Abstract: The safeguarding of resources is one of the principal subjects of 
haheutics studies. Among the solutions proposed to avert the disappearance 
of species, the setting in place of no take reserves is often mentioned. Most 
work on this subject, theoretical as well as applied, was undertaken in recent 
years. In this paper, we seek to compare two different models presented in 
existing literature by highlighting their underlying assumptions. Both models 
were derived from what is often referred to as the "model of Schaefer-Clark" 
(reference to the work of the last author on Mathematical Bioeconomics : 
Clark [7]). We show that various variations of this model lead to properties 
that can be very different. 

Keywords: dynamical system; calculus of variation; infinite horizon; ma- 
rine reserve; bioeconomic model 

1 Introduction 

It is now well known and largely accepted within the scientific community 
that the exploitation of halieutics resources has reached a critical thresh- 
old and solutions must be found urgently to conserve marine biodiversity - 
and indeed the existence of certain species. Beyond the application of quota 
policies, other strategies have been proposed such as the creation of marine 
reserves. The study of the role of reserves in fishery management has been 
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the subject of renewed interest in recent years [21 El El 1101 [HI [II]- The In- 
ternational Conference on the Economics of Marine Protected Areas (MPA) 
held in July 2000 in Vancouver was one of the starting points in the de- 
velopment of this new paradigm: the use of MPAs as an instrument in the 
management of fisheries. 

In economics literature, Sanchirico and Wilen [IB] seem to have been the 
first to suggest that MPAs could be beneficial not only from an ecological 
but also from an economic point of view. In their dynamic and spatial model 
of a Marine Reserve Creation, they analyze whether the transfer of biomass 
from the reserve to areas where catch is allowed could create economics profit 
from the MPA since its creation could help to improve a depleted biomass 
and increase catch outside of the reserve. They call this a double-payoff 
because in this case the MPA would increase both biomass and economic 
profits from the fishery. 

Both theoretical and applied aspects of the subject are well documented 
in the literature. From a theoretical point of view, numerous types of models 
have been proposed based on differential equations, mixing an inter temporal 
dynamic that corresponds to the population growth under consideration with 
a spatial distribution of this dynamic over diverse zones [9l [13 [H] . 

Sanchirico and Wilen construct a model where "the population struc- 
ture is characterized in a manner consistent with modern biological ideas 
that stress patchiness, heterogeneity and interconnections among and be- 
tween patches" [IT]. In this model, independent growth dynamics thus are 
associated with different patches. 

Another type of model has been used in the literature to account for an 
analogue structure [19] where a population develops different characteristics 
in sub-zones. In these models, the population of the entire zone under con- 
sideration follows a given dynamic evolution and the diverse sub-zones have 
dynamics such that by aggregating them together one may rediscover the 
global dynamic. 

The question that one then must ask is whether these two approaches 
may be used interchangeably. Few studies have focused on comparing these 
different model types. We would like to demonstrate that the choice between 
a patch model and a global model is not a neutral one, and highlight char- 
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acteristics of these two models that often are not specified in the hterature. 

The two models that we will compare are both Clark type (Schaefer, 
Gordon,...) [7] whose dynamic is logistic (Verhulst) and therefore widely 
used in halieutics dynamics. We use them to study the consequences of 
setting up a marine reserve from both an economic (inter temporal revenue) 
and a biologic point of view (population stocks). 

The first, the patch model, assumes a relative autonomy between reserve 
and non-reserve zones. This type of model is fairly widely used. The second 
model, known as the global model, assumes for its part a greater interaction 
between zones. 

The rest of the paper is organized as follows. In Section 3, we introduce 
and examine the two models we wish to study. We then compare results 
obtained, particularly using numerical simulation, in Section 4. Section 5 
concludes and is followed by a series of annexes that present the demonstra- 
tion of various results. 

2 Two variations on the Clark model 

We present two modelizations for a protected area in a given zone. These two 
modelizations derive from the well known fishery model studied for instance 
by Clark |^ among others. 

2.1 The first variation: patches model 
2.1.1 The model 

We consider a fish population that lives in a zone caracterised by a car- 
rying capacity K = 1. We assume that this zone splits in two sub- zones, 
with capacity respectively equal to a and 1 — a; we denote the stock of the 
corresponding sub-populations by Xi and X2- 

These two populations follow two independent evolutions laws and this 
is the reason why we use the "patches-concept" . These evolutions are given 
by: 

^(t) := x,{t) = F,{x,{t)) 
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The standard reference for the evolutions law is the logistic law: 

F,{x^) - r^x^{l-^) 
a 

F2{x2) := r2X2(l - , )■ 
1 — a 

But for our purpose it is enough to assume that the Fi are strictly con- 
cave, functions defined on [0, a] respectively on [0, 1 — a] and the Fj satisfy 
F,(0) = 0,Fi(a) = F2(l-a) = 0. 

We now assume that some exchange exists between these two patches 
and that this can be represented in terms of the density of these populations. 
More precisely we assume that the existence of some diffusion between these 
two patches can be captured by the following: 

1 — a a 

where A > represents a diffusion coefficient (3^, ^ being the density of the 
populations). The value of the diffusion coefficient depends on the location 
of the protected area. 



Prom now on, we decide that the first zone with capacity a is a protected 
area where no catch is allowed, whereas in the second zone fishing is allowed. 
The "normal" situation, i.e. the protected area acts like a source of biomass, 
corresponds to the case where the density inside the protected area is bigger 
than outside, i.e. 

X2 ^Xi 



1 — a ~ a 

The growth of the two sub-populations are governed respectively by the fol- 
lowing dynamics: 



x,{t) = Fi(xi(t)) + A(t)(^-^) 

Ht) = F2{x2{t)) - \{tr^ - ^-^) - h{t) ^' 
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where h(t) is the capture rate at time t. 

We note, from the positiveness of A and of the functions Fi{.), that if the 
sytem ([T]) possesses an equihbrium, it has to be necessarily normal 

As it is generally assumed the catch is proportional to the fishing effort 
E, and to the density of the population [1], therefore given by: 

hit) = qE{t)^ 
1 — a 

where q, the catchability coefficient, represents the fishing death rate when 
the density of the population is equal to one. We assume 

< E{t) <Em q>0 



The catch is sold on a market. In order to simplify we assume a constant 
price, p, over time and a constant cost, c, proportional to the effort. Therefore 
the revenue at t time is 

phit) - cE{t) = {pq^ - c)E{t) 
1 — a 

We then consider the discounted total revenue on an infinite horizon is 
given by 

J(E(.), A(.)) := r e''\pq^ - c)E{t) dt (2) 

where 5 > is an actualisation factor. 

We assume the existence of a manager whose goal is the maximisation 
of this total revenue. Moreover we assume that this manager can act on 
the fishing effort E and on some caracteristics of the reserve (closure, loca- 
tion) that are captured by A. Then the manager faces the following control 
problem: 



max JiEiXXi.)) 
E{.),H.) ^ ^ 

s.t. ([ID (3) 

< E{t) < Em 



Remark In many papers the states variables stand for the densities of 
the populations and not for the amount of the biomass. The link with the 
present model is obtained in setting: 

_Xi X2 
^1 — 5 ^2 



a 1 — a 

The two dynamical equations that give the evolution of the populations in 
the logistic case, then become: 

X^ = riXi(l-Xi) + ^(X2-Xi) 

X2 = r2X2{l-X2)-j^{X2-X,)-QEX2 

with Q = j^. 

This model with patches could be considered as a more general one than the 
P paper which corresponds to a = ^. 

2.1.2 Analysis of the solutions 

We will study the previous optimal control problem by the help of the calculus 
of variations theory . 

Prom the dynamic ([1]) , we deduce the expression of the effort in terms of the 
state variables: 

1 - a 

E{t) = -r(Fi(xi(t)) + F2(x2(t)) - xi{t) - Ht)) 

qx2{t) 

and then we obtain the new form of the objective. Thus the optimisation 
problem becomes: 



^ s,, c(l-a] 



max/ e (p )\Fi(xi) + F2(x2) - xi - X2] dt 

^GC- Jq qx2 

C being the set of admissible curves: 

C = {x{.) = (a;i(.), X2{.)) x,{.) G ^^^([O, oo[), x,(0) given , 

^(xi, X2) - qEuj^ <xi + X2< G{xi, X2)] ^ ' 

with G{xi,X2) := Fi[xi) + F2{x2) and BC^ stands for the bounded with 
bounded derivative functions defined on the interval [0, 00 [. 
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It is known that on BC^, the first order optimahty conditions given by 
the Euler-Lagrange equations, apply (see [3]). We suppose that x{.) stands 
for an interior solution and then x{.) has to satisfy: 

/^,(xi(t),X2(t)) - 4/i,(xi(t),X2(t)) +5/i,(xi(t),X2(t)) = 

at 

/(.,.) being the non actuahsed Lagrangian of the calculus of variations prob- 
lem, IxX-y) stands for the derivative with respect to Xi and IxX-, ■) stands for 
the derivative with respect to Xi.. 
The Euler-Lagrange equations becomes then 



*i = x,{^^-l){Fi{x,)-S) + F^{x^) + F,{x,) 
^2 = x,{^^-l){6-Fi{x^)). 

We first are interested by the non trivial equilibria , {xl,xl), of ([6]) i.e. such 
that Fi{x*) 7^ 0, i.e. xl 0,a and X2 7^ 0, 1 — a. It is easy to estabhsh that 
such equilibria have to satisfy (Appendix 1) 

F[{xl) = 6. 

We now assume that ri > 6. Then in the logistic case, from the strict con- 
cavity of Fi{.), we immediatly deduce the existence of a unique g]0, a/2[. 
Therefore we obtain the following result whose proof is postponed in Ap- 
pendix 1. 



Lemma 1 In the logistic case with ri > 5, there is a unique positive non 
trivial solution, {xl,X2), of the Euler-Lagrange equations This solution 
is caracterised by 

, a{ri - 6) 

and X2 given by 

c(l — aV cil — a) 1 — a 4ri 



7 



Clearly x* G]0,a[. It remains to show that G [0,1 — a]. This can 
be done straightforwordly but we prefer to use the following approach. We 
recall that if ([1]) possesses a non-zero equilibrium (x^, X2) then it is necessarily 
normal that is to say < — and we first observe that if this last condition 
holds then Xg < ^(1 — a) < 1 — a. 

We prove now that this condition holds under conditions on the ratio ^ and 

on the size of the different zones given by a. 

Let us introduce T(.) defined from the left hand member in ([7]) 

T{z) = z[2r29z^ - {9{r2 - S) + r2)z - S] 

where 6 := —. Then: 

c 

x*2 X _ g {ri-S){ri + S) 
1 - a 1-a 4ri 

As we assume ri > 6, then T(-^) > 0. From the graph of T(.) we deduce 
that the inequality 

o 1 



1 — a a 
holds if and only if 

T(-^)<T(^). 
I — a a 

Let us introduce when ri 7^ r2 

«.:=(J:i^ + ^)^. (8) 

ri — Ti — r2 

The following result is detailled in Appendix 2. 

Lemma 2 

/n i/ie logistic case, 

1. if ri> r2, for each 6 > Oq if ot satisfies 

a{n + 6 + '-1^(66'-^ - r,) - 26) < '-^{96'-!^ - r,) - 26 
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then the solution {xl,X2), of the Euler Lagrange equation ^ given in 
Lemma 1, is normal i.e. satisfies: 



— a ~ a 



2. if Ti < r2, (x^jXg) is never normal i.e. 



a a 



Then it remains to prove that ( candidate to the optimisation 

problem ([3]), thus that the constraints are satisfied. 

If the conditions in the Lemma 1 and Lemma 2.1 hold, then Fi{x\) > 

and — - — < — . Therefore to this corresponds a unique coefficient of 

1 — a a* 
diffusion A* > 0. 

From (II]) we deduce that 

E* = ^{F,{xl) + F,{x;)) > 0. (10) 

Moreover the expression of the total revenue is given by 

r = {pqj^-c)^. (11) 
1 — a 

This revenue is positive if the fishery profit is positive at this equilibrium, 
that is to say if 

X2 c 1 

1 — a pq 9 
But this inequality holds because we have 

and due to the fact that 6 > 6o > ^, this last term is nonpositive . 
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Proposition 1 

In the logistic case if ri > 5 the problem ([^j possesses at most a non trivial 
and positive optimal stationary solution caracterised by 

^ a{ri-5) 

X, = 

^ 2ri 

and X2 given by The corresponding effort, diffusion coefficient and total 
revenue are given respectively by (TB^) . (Cp, ( fii]) . 
Remarks 

1) When Oi = \ the value of 6q coincides with the value Pm given in [1]. 

2) From the expression of 6*0 in ([8]), we observe that ri can't be closed to r2- 
If this is not the case, then the value of the dimensionless ratio 9 has to be 
very high. But this can be unrealistic because the value of 9 is given by the 
economic environment. 

3) From ([8] ) with a given value for 9 we can precise a bound for r2 expressed 
in terms of ri, 5: 



9 



2ri 



To < ri . 

2 - Q^r^ 

2.2 The second variation: the sphtting of a unique 
zone 

In this second model we start with an unique zone with capacity K that 
we normalise to one, K = 1. Let us assume that the population follows a 
standard evolution law: 

z{t) = (f){z{t)) (for instance = rz{t){l — z)) 

being a concave function defined on [0, 1], with 0(0) = 0(1) = and 
0'(1) < 0. We assume that this zone splits first in a part that is a reserved 
area where no fishing is allowed and a complementary part that is open to 
harvest. We assume that these two parts have respectively a and 1 — a as a 
carying capacity. 

The main difference with the previous model is that the two populations. 
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which stocks are respectively Xi and X2 follow the evolution laws given by: 

Xi{t)=F{x^{t),Xi{t)+X2{t)) 
X2{t) = F{x2{t),Xi{t) + X2{t)) 

where F{., .) satisfies the standard assumption of regularity with F{x, z) — 
if and only if x — or z — 1 and where F(.) and 0(.) satisfy 

F{xi, xi + X2) + F{x2, xi + X2) = (j){xi + X2). 

For instance F can be a logistic function 

F{x, z) = rx{l — z) = rx{l — {xi + X2)). 

As in the previous model, there is some diffusion between the two zones which 
can be represented by: 

1 — a a 

Then the two populations evolve following the dynamics: 

x,{t) = F{x,{t),x,{t) + X2{t)) + - ^) 

1 — a a 

X2{t) = F{x2{t),X,{t) + X2{t)) - - 

1 — a a 

We want to stress on the fact that this new model is consistant in the sense 
that the sum of the two dynamics is exactly the evolution law of the total 
population. 

Now taking into account the catch in the zone where fishing is allowed, 
we derive the final dynamics of the populations 



x,it) = F(xi(t),,xi(t) + ,X2(t)) + A(^-2^) 

X2{t) = F{x2it),x^{t)+X2{t))-XC^-^-^)-qE^^ ^ ^ 

where E e [0, Em] stands for the fishing effort and g > is the catchabihty 
coefficient. 

In order to compare with the previous model, we assume that a manager 
has as an objective to maximise the actualised total revenue as presented 
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before. To do so, he has to act on two controls, the fishing effort and the 
location of the reserve area given by A. Therefore the manager faces the 
following program of optimisation 

max r e'^Hpq^ - c)Eit) dt 

E{.), A(.) ^ ^ ^ 

s.c. m (13) 

< E{t) < Em 



2.2.1 Analysis of the solutions 

From the dynamic equations f[T^ we can derive the expression of the effort 
1 - a 

E{t) = 7-(0(a;i(t) + X2{t)) - xiit) - X2{t)) (14) 

qx2{t) 

and then we obtain the equivalent problem to (fT3|) as a calculus of variations 
problem: 

max / e [p j\(p{Xi + X2) — Xi — X2\ at 

xec Jq qx2 

where C stands for the set of feasible curves defined by: 
C = {x{.) = (xi(.), X2{.)) Xi{.) e BC\[0, oo[), Xi{0) given , 

0(Xi + X2) - qEufl^ <Xi+X2< 0(Xi + X2)}. 

We know that in this framework a necessary optimality condition for an 
interior solution x{.) is given by the Euler-Lagrange equations that are: 

ii = X2(^-l)(0'(a;i+a;2)-5) + 0(xi + X2) 

±2 = X2(^-l)(5-0'(Xi+X2)). 

For now, we will stick to the logistic case. The Euler-Lagrange equations are 
then: 

^1 = ^2(^g!% - 1)('^ - 2r(xi + X2) -S)+ r{xi + X2){1 - (xi + X2)) 

X2 = X2(^-l)(5-r + 2r(xi+X2)). ^^^^ 
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In order to derive the non trivial positive equilibria, denoted by (x*,X2), 
we first consider the second equation in f|T5|) with the condition 



This implies result 

= r(xi + X2)(l - (xi + X2)) = — . 

A contradiction if r > 6. In the case where r = 6, then we obtain the trivial 
solution xl = X2 = 0. 

Then we deduce that an equilibrium has to necessarily satisfy 

_ c{l-a) 
pq 

With the help of first equation in (fT5|) we find that either + = or 

Xl+ X2 = 1. 

Finding a non trivial equilibrium implies to exclude the first condition. There- 
fore we have proved the following result 



Lemma 3 In the logistic case there is a unique non trivial and positive 
solution for the Euler- Lagrange solutions ( [73]) given by: 

f * *^ c{l-a) c(l-a) 

Xi,X2 = 1 , . 

pq pq 



In order to examine whether this candidate solution of the problem (I13p 
can be optimal or not, we have to derive the corresponding effort and diffusion 
coefficient. From the expression of the effort dHl), we obtain 0(x* + X2) = 
0(1) = that is to say 

E* = 0. 

We also deduce that 

A* = 
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except if ^ = 1. Finally at this equilibrium the intertemporal revenue is null 
too. We have established the following proposition 



Proposition 2 

In the logistic case the problem /[T^) possesses at most a non trivial and 
positive stationary solution given by 



pq pq 

The corresponding effort, diffusion coefficient and total revenue are null. 
Remarks 

1) It is easy to obtain that this equilibrium is normal, i.e. 

— > 



a 1 — a 
if the fishery is profitable, that is to say if 

pq — c> 0. 



2) An adaptation of the model given in Gomez et al. [6] to our case of a 
no take zone is: 

= ar{xi + X2){1 - {xi + X2)) + - , . 

X2 = {l-a)r{x, + X2){l-{x,+X2))-\{^-^)-q2E2^. ^'^^ 

We can get the same results as those given earlier. In Gomez et al. [6], fishing 
is allowed in the so called artisanal zone (corresponding to the protected area 
in our case) and the objective to be maximised is somewhat different (it takes 
into account the revenues of the artisanal and industrial fisheries). Here also 
it has been proved that a unique solution exists but with a non null effort 
and a non null revenue. 
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3 Comparison, Numerical application 



In this section we underline the differences between the results we obtained 
in the previous sections for both the patches case and the global model. 
Prom their expressions given in the Propositions 1 and Proposition 2, we can 
make the following remarks for the equilibria {xl,X2): 

• In the model with patches the first component doesn't depend ex- 
plicitely on the ratio ^, whereas it does in the global model. 

• In the global model the second component is given by X2 — 
whilst the patches model doesn't possess any equilibrium with such a 
component as the second, cf. Appendix 1. 

Thus the expressions of the equilibria are different in the two models. 

Now we estabhshed that the optimal effort and the corresponding total 
revenu at (a;*, X2) was null for the global model. This doesn't seem to be the 
case for the patches model, we will later show with simulations that optimal 
effort and total revenue are not significantly close to zero. 

In order to continue the comparison, let's arbitrarily fix the parameters 
p, q, c, a, S. Thus the models depend only on the instantaneous growth rates 
ri,r2 and r respectively. 

If we let ri — r2 — r, in Lemma 2 we established that the equilibrium {xl, X2) 
was never normal in the patches case, while it is always normal for the global 
model (Remark 1 in section 3.2). 

Now to compare our models with ri and r2 only near r, we noted in the 
Remark 2 of section 3.1 that this situation wasn't a realistic one. 

Then comparing these two models is not an easy task. Our first conclusion 
is: the role of the instantaneous growth rates of the biomasses are crucial to 
choose such or such model. An assumption that is not underlined in general. 

Now let's come back to the comparison of the optimal efforts and revenues 
by using simulations. The main issue is to determine significant growth rates 
that are not equal. 

But this choice shouldn't depend on our particular models with preserving 
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areas. It should be the same for a wide class of models. For instance for 
models that correspond to a situation where fishing is allowed in the the two 
areas ([13]). We consider thus: 



Xl = Fi(xi) + A(^ 

X2 = F2{X2)-X{^ 



a 



a 



) — qExi 
) - qEx2 



(17) 



and 



Xl = F{xi,xi + X2) + \{ 
±2 = F{x2,xi + X2) - X{ 



X2 
1-a 

1-a 



Xl 

a 

Xl 



a 



) — qExi 
) - qEx2 



(18) 



with the same assumptions as before. In order to compare numerically ( ITTll 
and (|T8|) we will face the same issue to determine significant growth rates. 
We propose to use this new situation in order to fix values for ri,r2,r. The 
new problem we consider now is to maximize the same objective as before 



but with ( IT71) and ( fT8l) . We observe that (|T8l) corresponds to the classic 
Clark model, it is enough to let z = Xi + X2 to obtain that the dynamic is 
i = (p^z) — qEz. 

We can assume that a manager has no reason to use one model rather 
than another. Then the two models can be considered as equivalent in the 
sense that they provide the same optimal effort. 

Hence, we propose the following procedure to determine a system of growth 
rates: Let's set an arbitrary choice of values for ri and r2- From the first 
order optimality conditions, given here by the Pontryagin principle, we can 
derive the optimal value of the corresponding effort for problem (fT71) . We 
hand-over this value in the first order optimality conditions of the second 
problem ( fTSl) from which we derive the value of the growth rate r. 
To follow this procedure we set: 




a = .5, 5 = .05, c = .15, g = 2., A = 20 



and we obtained for 



n = 0.4, rs = 0.05 
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that 



E = .0566 and r = 0.28739. 



Let's now go back to our models with protected areas from where we take 
the previous values for the parameters and where we set ri = .4, r2 = .05 
and r = 0.28739 . 

Then for the model with patches we found that: 

• the optimal effort is E* = .0457, and the biomass values are respectively 
xl = .21875, = .0302 

and for the global model 

• E* = Oetxl = .875, = .125 . 

We observe that the optimal efforts corresponding to the patches case, 
E* = .0457, and the Clark model ( ITH]) . E = 0.0566, have similar sizes. We 
know that the optimal value of the effort in this last model can't be considered 
as null. Therefore we can deduce that in the first model with patches the 
effort is not null. 

Then the two models have different qualitative behaviour:total revenues and 
optimal effort are totally different. 

4 Conclusion 

In this article, we have shown that different models have been proposed 
and used in the literature for the same MPA problematic. We focused on 
demonstrating the importance of the hypotheses underlying two types of 
models - the patch model and sub-zone model -particularly the crucial role 
played by the growth functions (rate and form), and on studying the different 
results produced by them. 

A manager who wishes to study the role of an MPA in a given zone 
must first know if the entire zone is artificially divided or if it can be broken 
down into patches (entities with their own dynamics). Without taking this 
precaution, and in obtaining the very different results that we have seen, the 
manager risks taking erroneous decisions. 
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The two preceding models of resource dynamics are adapted to the case 
where control instruments are independent of the size of the no take reserve. 
If the manager must take size into account in his decisions, the modelling 
of the dynamic has to be changed. For instance, it is necessary to consider 
a depending of a diffusion coefficient. A justification is given in Appendix 
3. In this Appendix, we also underline that this coefficient can be given by 
formula 

\{ol) = Aoa;(l — a) 
which is the expression considered by Boncoeur (cf. 

5 Appendix 
5.1 Appendix 1 

We determine the non trivial equilibria of the Euler-Lagrange equations ([6]) 

^1 = X2{^^ - im{x2) - 6) + F,{x,) + F,ix,) 

If we assume that = p""'' , from the first equation: 

F,{xi) + F^ixl) = 0. (19) 

From the assumption of the non triviality of the equilibria, we have that 
c 7^ 0, a 7^ 1, c 7^ and then -^2(3^2) > 0- Therefore we can't find any 
Xi G [0, a] such that ( TT9l) holds. Thus there is no non trivial equilibrium 



with X 



2 



c(l-a) 
pq 



Therefore in order for a non trivial equilibrium to exist it is necessary 
that 

FiixD) = 5. 
In the logistic case it is easy to compute that 

a(ri — 6) 



x^ 



2ri 
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Reporting this value in we get the following equation for 

c[l — a) 
which is ([7]) in the logistic case 

2r2pq 2 ( ( ex N ,--1 ('^i - 5) (ri + 5) 



X2 



71 \^^2 - (^1 ^(^"2 + )X2 -6] = a- 

c(l — a)^ c{l — a) I — a 4ri 

The graphes of the functions defined by the left and right hand members 
are curves that crosse in a single if ri > 5. But this last inequality holds 
because from our assumption we have 

p. *^ (n-5)(ri+^. n 
^i(^i) = « > 0- 

This ends the proof of Lemma 1. 
5.2 Appendix 2 

In order to find conditions for the inequality < 77 to be true, we know 
that it is equivalent to consider the inequality < T(^). This last 

inequality becomes 

a (ri-3)(ri+3) ^ 
1— a 4ri — 

that is equivalent to 

a(ri + 5) < (1 - a)C-^^{e5'^^—^ - r^) - 25). 

ri ri 

1) If ri > r2 the right hand member has to be positive, this implies the 
following condition on 9: 



ri — 6 5 ri — r2 
Now if this condition on 6 holds, from the previous inequality we should 
deduce (E]). 

2) If ri < r2, the right hand member is always negative and therefore 
(^X-^ , can 't be normal. 

This ends the proof of Lemma 2. 
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5.3 Appendix 3 



We consider the case where the manager has the size of the preserved area 
as control. We will first prove that the diffusion coefficient has to depend on 
this size. 

We start with the dynamics and the objective given in the second variation 
(§ 3.2). We suppose that the manager has to maximise his objective by using 
the fishing effort E and the size of the preserved area that is captured by a. 
We always denote by z the stock of the total population and the two sub- 
populations stocks by xi and X2 respectively. Then the densities in the two 
regions are di = — and d2 = -p-. 

When a = we can only find a single zone and thus Xi = and X2 = z. In 
this case, it is natural to set for the densities: cii = and d2 = z respectively. 
Now if a = 1, it is natural to set: di = z d2 = 0. 

As we have done before, we assume that some diffusion exists between the 
two zones and that it is proportional to the difference of the two densities. 
Therefore in order to respect our previous remark, we have to set 

A(«)(t^--) (20) 
1 — a a 

where the diffusion coefficient depends on a. Indeed if a = 0, from di = 
we deduce that A(0)^^ = X{0)z = because in this case we can only find a 
single zone, and thus A(0) = 0. From a similar argument, we deduce that for 
a = 1, we have A(l) = 0. Now for a g]0, 1[ the coefficient X{a) is certainly 
not null. 

For instance we can let 

A(a) = Aoa(l -a) (21) 

and in this case the diffusion is modelised by 

Aoa(l-a)(-^-^). (22) 
L — a a 

This expression is the one proposed by Boncoeur in [1]. 
Then the problem of the manager becomes in this setting 

max r e-^\pqf^ - c)E{t) dt 
s.c. ([T2]) 
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where the diffusion coefficient in fll2l) is given by fl21l) . 

The solution of this problem is straightforward using the Pontryagin maxi- 
mum principle. We won't mention it in this paper. 
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